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Abstract 

The target of this work is to investigate the physical nature of polar jets in the solar corona 
and their possible contribution to coronal heating and solar wind flow based on the analysis 
of X-ray images acquired by the Hinode XRT telescope. We estimate the different forms of 
energy associated with many of these small-scale eruptions, in particular the kinetic energy 
and enthalpy. Two Hinode XRT campaign datasets focusing on the two polar coronal holes 
were selected to analyze the physical properties of coronal jets; the analyzed data were 
acquired using a series of three XRT filters. Typical kinematical properties (e.g., length, 
thickness, lifetime, ejection rate, and velocity) of 18 jets are evaluated from the observed 
sequences, thus providing information on their possible contribution to the fast solar wind 
flux escaping from coronal holes. Electron temperatures and densities of polar-jet plasmas 
are also estimated using ratios of the intensities observed in different filters. We find that the 
largest amount of energy eventually provided to the corona is thermal. The energy due to 
waves may also be significant, but its value is comparatively uncertain. The kinetic energy 
is lower than thermal energy, while other forms of energy are comparatively low. Lesser 
and fainter events seem to be hotter, thus the total contribution by polar jets to the coronal 
heating could have been underestimated so far. The kinetic energy flux is usually around 
three times smaller than the enthalpy counterpart, implying that this energy is converted 
into plasma heating more than in plasma acceleration. This result suggests that the majority 
of polar jets are most likely not escaping from the Sun and that only cooler ejections could 
possibly have enough kinetic energy to contribute to the total solar wind flow. 
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1 Introduction 


Coronal jets are small eollimated ejeetions of plasma observed in ultravio let/x- 
ray imagers and in white-light coronagraph images [see St. Cyr et^D. 1997 1. The 
kinematies of polar jets were first studied by Wang et aP [ 1998 1 during solar mini- 


[LASCO; 

Brueckner et ah. 

1995 

] and the Extreme ultraviolet Imaging Telescope 

fEIT: Delaboudiniere et ah. 

I 995 I 1 onboard the SOlar and Heliospheric Observa- 

tory [SOHO; 

Domingo et al. 

. 19951. The authors analyzed 27 correlated events and 


derived outflow speeds of ~ 250km ■ s~^. More reeentlv. ICirtain et al.l OOOVll used 
Hinode IIKosugi et al.l. 1200711 observations of polar eoronal holes to reveal that x-ray 
jets have two distinet veloeities: one near the Alfven speed (600 - 800 km ■ 5“^) and 
another near the sound speed (~ 200km ■ 5'“^). This indicated that the jets may con¬ 
tribute to the high-speed solar wind. These observations also demonstrated that jets 
are transient phenomena that occur at much higher rates than large-scale events, 
such as flares and coronal mass ejections. 



board Hinode showed that around 60 x-ray jets per day occur on average inside 
the polar coronal holes. The apparent outward velocity, the height, the width, and 
the lifetime of the jets were measured, and a statist ical study of jet transverse mo¬ 


tion s and backflows was also pre sented. Recently, IChandrashekhar et al.l rt2014all 


and IChandrashekhar et al.l 1120 14bn studied the dynamics of two jets observed with 


the Hinode/XRT and Hinode/EIS and confirmed the drift motions that appear in 
the upward movement of the event, suggesting that the jet material is undergoing 
cooling and falling back. 

Correlations between jets and ot her signs of solar activity have been studied, for 
instance, by Raouafi et al. I 2008ll . who focused on the temporal evolution of and 
relationships between polar coronal jets and polar plumes, a nd showed that ~ 90% 


of the observed jet events are associated with polar plumes. iFilippov et al.l [l2009ll 
discuss the formation of jets and proposed scenarios that explain the main features 
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of the events: the relationship with the expeeted surfaee magnetism, the rapid and 
sudden radial motion, and possibly the heating, based on the assumption that the 
jet occurs above a null point of the co ronal magnet i c field . Using the Solar Terres¬ 
trial Rel ations Observatory [STEREO: iKaiser et all 1200811 twin inner coronagraphs 
rCORl: Thompson et all 2003 1 Nistico et al. I 2009ll studied various aspects of jets, 
including their correlation with un derlying small scale chromospheric bright points. 
More recently, Pucci et ^ I 2012ll have presented an analysis of the correlation of 
X-ray bright points and jets observed by Hinode/XRT and concluded that most of 
the jets occurred in close temporal association with the brig htness maxirnum in 
brigh t points. A dichotomy of polar X-ray jets was proposed bv iMoore et al.l 112010 . 


201311 . They observed that about half of their sampled events fit the standard re¬ 


connection picture for coronal jets, and the rest correspond to another type. The 
non-standard jets (which they call "blowout jets") are described as having a jet- 
base magnetic arch that undergoes blowout eruptions, similar to, but smaller than 
those that produce major coronal mass ejections. The authors also propose possi¬ 
ble correlations between jets and Type-II spicules and macrospicules and conclude 
that the combined energy output could contribute significantly to the heating of the 
corona. 


The ph ysical parameters of jet plasma have been deri ved, for instance, by 


[ 2 OOOI 1 with SXT data on Yohkoh. More recently, Nistico et al. [ 2011 1 have de- 


Shimojo and Shibata 


scribed the typical physical characteristics of coronal jets observed by the SECCHI 
instruments of STEREO spacecraft obtaining a temperature determination for the 
jet plasma. Their results show that jets are characterized by e lectron temperatures 
rangin g between 0.8 MK and 1.3 MK, similar to results by lYoung and Muglach 


[ 2014 1. who studied the link between a blowout jet and its base brightpoint. lPucci et al. 


[1201 311 analyzed the difference between a standard and a blowout jet seen from 
Hinode-XRT. They aimed at inferring dif ferences in physica l parameters corre- 
sponding to the two categories proposed bv iMoore et al.l KOlOll . lPucci et al.l [l2013ll 
conclude that their blowout jet was hotter, had a higher outflow speed, and was 
rooted in a stronger magnetic field region compared to their standard jet. 


In our previous work [ Paraschiv et ah . 2010ll . we identified and studied white- 
light jets observed by the SECCHI-CORl white light coronagraph. We identi¬ 
fied more than 10 000 white-light jets spread across the entire solar limb at the 
minimum of solar activity between 2007 and 2008. The identified jets originated 
in regions inside coronal holes but also in regions of the quiet Sun. A subsam- 
ble of the events were considered for correlation with bright points on the solar 
disk. The association was done based on the changes (morphological, brighten¬ 
ing, disappeari ng, etc.) of the brigh t points as observed in Extreme Ultraviolet 


Imager [EUVI; IWuelser et al.L 1200411 195A images at the time when the jet was 


observed in CORl. In Paraschiv et al.l [ 2010ll . we suggested that the number of jets 
with high outflow speeds is large enough to contribute significantly to the solar 
wind, and we concluded that these results may point to the mechanism responsi¬ 
ble for the expansion of the jet, namely a pressure-driven expansion of the plasma. 
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We later expanded the study rtParasehiv et all 1201211 in order to diseuss the pos¬ 
sible contribution of coronal jets to the solar wind flux, by estimating the ejected 
particle mass flux of a selected subsample of the above-mentioned coronal jets. 
In the present work we want 
to extend our previous analy¬ 
ses to X-ray data and to de¬ 
rive plasma physical parame¬ 
ters on a statistically signifi¬ 
cant sample of jets. The data 
that were used, along with the 
physical phenomena and mod¬ 
els and assumptions involved in 
the analysis, are described in 
Section 2. The third section de¬ 
picts the methods used, compu¬ 
tations, and results obtained: the 
jet’s geometrical and ki nematic 
properties, jet temperature, and 
density determinations are de¬ 
scribed. An estimate of the en¬ 
ergy budget for coronal jet erup¬ 
tions is also presented. In the 
fourth Section, discussions, re¬ 
sults and conclusions are presented. Comparisons to previous studies are given and 
conclusions about possible contribution to coronal heating are made. 



response, 
filters we 


Fig. 1: Total XRT temperature 
given for the X-ray focal-plane 
used. Each curve plots the combination of 
the total instrument response as a func¬ 
tion of wavelength, assuming a coronal 
plasma emission model (ATOMDB/APEC). 
The full filter response curves are available in 


Narukage et al.l ll2014ll . 


2 Data selection & calibration 


All the data analyzed in this work were acquired by the XRT instrument onboard 
the Hinode spacecraft. The XRT is a modified Wolter I telescope that uses grazing 
incidence optics to image the solar corona’s emissions (T>~1MK) with an angular 
resolution of 1 arcsec per pixel at the CCD. The telescope has a maximum imaging 
field of view of 34 arcminutes. It uses several filters t uned to differ i ng wa velengths, 
allowing sampling of different regions of the corona. iGolub et al.l 11200711 gives full 
details for the instrument. 


Figure [H shows the response functions, FfuiTg), for the filters used here as func¬ 
tions of temperature. XRT operates on a mission schedule and changes its pointing 
accordingly. For the purpose of this analysis, we had to search for suitable sets of 
data. Instrument pointing in polar coronal holes was required for a good (via a better 
signal-to-noise ratio) differentiation between the jets and the coronal background. 
Full resolution of 512 x 512 pixels with a spatial resolution of 1.02 arcsec/px and 
adequate image cadence (temporal resolution) of no longer than ~ 2 min were used 


4 












as mandatory data conditions for this study. The data also needed to be available in 
multiple filter eombinations in order to apply the filter ratio teehnique used to de¬ 
termine parameters of the jet plasma, sueh as the temperature, density, and energy 
budget of events. 


We seleeted the following datasets: 


• 1. Dataset aequired on 6 
July 2008, eontaining approx¬ 


-800 

imately four hours (10:36 to 
14:40) of eontinuous observa¬ 


-900 

tions of the north polar eoro- 

tn 


nal hole. Data were avail¬ 

< 

-1000 

able at full resolution, in three 
filters (Al_Poly, Al_Mesh, 


-1100 

C_Poly) with ~ 1.5 minute 
eadenee between eaeh three- 


-1200 

filter exposure set. 





-100 0 
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200 


Fig. 2: Example of XRT data field of view for 
the observations aequired on 26 Jan 2009. 


• 2. Dataset aequired on 26 
January 2009, eontaining ap¬ 
proximately five hours (07:02 
to 11:54) of eontinuous ob¬ 
servations of the south polar 

eoronal hole. Data were available at full resolution, in three filters (Al_Poly, 
Al_Mesh, Ti_Poly) with approximately a two-minute eadenee between eaeh 
three-filter exposure set. 

An example of a typieal XRT instrument field of view is shown in Figure |2l Data 
were retrieved online through the arehive seareh page of the Hinode Seienee Data 
Center (SDC) Europe (http://sde.uio.no/seareh/). All the data were analyzed using 
the Solarsoft (SSW) paekage, written in IDE language. The xrt_prep.pro proee- 
dure that was used to ealibrate raw XRT data from Eevels 0 to 1 mainly performs 
the following aetions: 1) eorreetion for undersaturated and oversaturated pixels, 2) 
removal of spikes that oeeur due to the radiation belt and/or eosmie rays, 3) eorree- 
tion of images for eontaminating material that has aeeumulated on the instrument, 
4) eo-alignement of the images by applying an orbital drift eorreetion to reduee 
the effeets of satellite jitter, 5) normalization of the output image in DN s“\ and 
6 ) subtraetion of instrumental dark noise using dark frames aequired elose in time 
to eaeh of the two datasets. Correetion for eontamination is very important; sinee 
Hinode’s launeh, eontaminating material has aeeumulated on the XRT CCD and 
foeal plane filters, eausing a deerease in sensitivity. The effeet of eontamination 
is wavelength-dependent, and the thiekness of the eontaminant deposited on the 
foeal plane filters is diffe rent, so that eontamiri ation affeets the response of eaeh 
XRT filter differently. See lNarukage et al.l 11201 ih for further details on analysis and 
modeling of eontamination effeets on XRT data and responses. 
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3 Data analysis 


3.1 Geometrical G kinematical parameters 


A total of 18 jet events were identified and seleeted by visual inspeetion of the 
approximately nine hours of observations that we analyzed. Image sequenees were 
employed to derive many different geometrieal and kinematical properties of these 
events. For each event, we first selected the frame with the jet at its maximum 
elongation prior to when it started to dissipate. A rotation was then applied so that 
the jet was vertical in our analysis window, and the region containing that vertical 
jet was preselected in a rectangular box. Then each event was isolated based on its 
intensity profile along the direction perpendicular to the jet’s motion; the resulting 
jet widths ranged from 3 to 12 pixels. 

Jet outflow speeds were estimated by summing each rectangular area over the 
direction perpendicular to the jet (jet width), providing the jet intensity distribu¬ 
tion along its axis, and then by tracking for each event the position of the leading 
edge, identified as the location where the intensity dropped below a fixed thresh¬ 
old, in successive frames. An example of a sequence acquired during one of the 
jets analyzed here is shown in Figure [3l Because all the jets we selected are ob¬ 
served off-limb, the plasma mainly propagates over the plane of the sky, so that 
the possible projection effects due to jet inclinations are minimized in the velocity 
measurement. For each individual event, data are available as independent series for 
each of the three filters, thus three outflow speed values were obtained. The speeds 
averaged over all the events are VAi_poiy = 157 + 31 km s~^,VAi_Mesh = 163 + 25 
km s“‘, and Vc_PoiyandTi_Poiy = 156 + 25 km s“‘; different outflow speeds for all 18 
events are provided in the jet parameter table (see Appendix). Computed errors are 
standard 1 cr deviations. More significant errors are due to uncertainties in the exact 
identification of the location of the jets’ leading edges, and the estimated resulting 
uncertainties in the outflow speed values are in the range of about 20%-25%. By 
summing each rectangular area over the direction parallel to the jet axis, we get a 
profile of the average jet intensity across its width: a 2cr dispersion from the loca¬ 
tion of the brightness maximum was selected as being the value of the typical jet 
width d. Values of the widths of the jets are also provided in the jet parameter table 
(see Appendix); the average widths of all of the jets is J - 8 ± 1 pixel (5900 + 750 
km). 


3.2 Plasma temperatures 


Using XRT data acquired with multifilter sequences, we can measure the electron 
temperature of the jet plasma via the filter-ratio technique. As mentioned, filters on 
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Fig. 3. Example of jet eruption in a sequence of base difference images (prior to rotating 
the frame for analysis; see text) showing the rising phase, the main phase, and the decay 
phase of a jet event. The outflow speed can be computed by time tracking the leading edge 
of the eruption. This sequence was acquired between 11:06 UT and 11:26 UT on 26 Jan 
2009 with the Al_Poly filter. 


the XRT telescope are characterized by a response function Ffu that depends on the 
temperature usually given in the units of [DN ■ 5“^ ■ pixeF^] for a unit column emis¬ 
sion measure (CEM) measured in [cm“^]. The measured intensity for a particular 
filter Ifii can be written as an integration along the line of sight (LOS) as 


Ifu= f Ffn(Te)nldl [DN ■ s-^ pixeF^ (1) 

Jlos 

The above equation shows that the measured intensities depend on the filter re¬ 
sponse, which is in turn a function of the plasma electron temperature (Tg). The 
filter response function also depends upon the abundance of the elements whose 
lines contribute to the emission measured by the instrument. 
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Fig. 4. Ratios computed with the filter response functions. 


The emitting plasma is eertainly not isothermal along the LOS. Nevertheless, be- 
eause the eleetron temperature distribution along the line of sight is unknown, some 
assumptions are needed to apply the filter ratio teehnique. In partieular, if the emis¬ 
sion due to the baekground eorona is subtraeted so that the emission due to the jet 
plasma alone is isolated, and if we ass ume the jet is isothermal along the LOS [see 
iMoreno-Insertis and Galsgaard . 2013 1. the above equation ean be rewritten as 


If a = FfuiT,) jnldl [DN • 5"' • pixel '^]. (2) 

The measured intensity also depends on the plasma eolumn emission measure, 
CEM = nl ■ I [em“^], where ne [em“^] is the eleetron number density. For this, 
we eonsidered the jets to have a eylindrieal geometry and thus the depth / is as¬ 
sumed to be equal to the projeeted width d on the plane of the sky of the studied 
events. 


Under our isothermal assumption the eolumn emission measure is the same for 
eaeh filter, and so the ratio R between the intensities measured in the same jet with 
two different filters (after subtraetion of eoronal baekground) solely depends on the 
eleetron temperature: 


IdTe) _ FdTe) _ DNdh _ ^ -^-i/ 7)iViAi \ 
h(Te) F2(T,) DN2/t2 " \DN2/t2l 


( 3 ) 


Equation [3] deseribes the direet dependenee between the emitting plasma eleetron 
temperature and the DNs deteeted for one jet with two different filters and different 
exposure times t, . Given the observed ratio and the theoretieal one, it is thus possible 
to derive the plasma temperature for eaeh pixel: different ratios employed for this 
analysis are shown in Figure |4] as funetions of temperature. As ean be seen in this 
figure, for some very limited ranges of filter ratios the temperature is not uniquely 
determined, beeause two different temperatures are assoeiated with the same single 
value of filter ratio. We handled these few double-valued situations by seleeting the 
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temperatures eloser to the average value obtained from the other filters, whieh is in 
general the lowest possible temperature. 

In summary, the following steps were followed in order to derive the average jet 
temperature. An intensity profile along the jet was ereated at the moment of its max¬ 
imum extent, together with a eorresponding intensity profile for the solar eorona 
obtained at the same spatial loeation and the same range of altitudes with frames 
aequired just before the jet’s start time. The latter profile was fitted with a seeond- 
order polynomial and removed from the intensity profile along the jet, thereby re¬ 
moving the estimated baekground and foreground eoronal emissions from it with¬ 
out inereasing the noise. From the ratio between the intensities observed with dif¬ 
ferent filters, the temperature was then ealeulated, using the eurves of Figure IH for 
eaeh ten-pixel interval along the jet. The resulting temperatures usually had roughly 
eonstant values along the length of the jet; as a result, by doing a global average 
over the entire jet, we get the average temperature of eaeh jet (Figure [5al). 


Only three filter ratios were eomputed for eaeh event owing to the different avail¬ 
able filter eombinations of the two used datasets. Some ratios give more reliable 
results than others: The derived temperatures are in general quite similar for the 
Al_Poly/Al_Mesh, Al_Mesh/Ti_Poly, and the Al_Mesh/C_Poly ratios, while the 
Al_Poly/C_Poly and Al_Poly/Ti_Poly ratio gives slightly different results. This is 
partly beeause some filters are more affeeted by eontamination than others, and also 
beeause the ratio of the response funetions of the Al_Poly/C_poly and Al_Poly/Ti_Poly 
filters is too fiat in that range of temperatures (see Figure H]). Considering that there 
were substantially larger variations in temperatures derived with the Al_Mesh/C_Poly 
or Al_Mesh/Ti_Poly ratios eompared to the Al_Poly/Al_Mesh ratio, we eoneluded 
that the latter gave slightly better results. This ratio is then the only one we used in 
further measurer nents. The estimated u neertainties eomputed using the error esti¬ 
mation given by iNarukage et al.l 11201 ill are on the range 10% - 15%. The average 
jet temperatures are 1.8 + 0.2 • 10^ K, temperatures measured for single events are 
shown in Figure [52 these results will be diseussed in the next seetion. 


3.3 Plasma densities 


Plasma temperatures obtained with the filter ratio teehnique were used in turn to 
derive an estimate of the plasma density. Values for the densities of both jet and 
surrounding eoronal plasmas are eomputed at the moment of maximum brightness 
during the jet evolution. In this work we define the jet density Uejet as the quantity 
to be added to the density of the surrounding eorona ne cor in order to get the density 
measured from the observed X-ray emission, ne_obs = ngjor + nejet- The observed 
jet intensity in a given filter (under the hypothesis of nearly eonstant temperature 
alon g the LOS) after s ubtraetion of the eoronal baekground is given by Eq. [2] [see 
also lPueei et al.l.l2013ll . As previously mentioned, the assumption that the jet has a 
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Fig. 5. Jet parameters derived using the three filter ratio combinations. Left panel (a): Jet 
temperatures (j-axis) derived for the 18 jets studied here (x-axis) using different filter 
ratios. Right panel (b): Jet densities (jy-axis) derived, using Eq. for all events (x-axis) 
using different filter ratios. Different filter pairs are depicted using different symbols and 
colors (as shown in the legend). Computed values usually correspond closely between the 
Al_Poly/Al_Mesh, Al_Mesh/C_Poly, and Al_Mesh/Ti_Poly ratios, while values derived 
from the Al_Poly/C_Poly or Al_Poly/Ti_Poly ratio are somewhat different. 


simple cylindrical form was adopted, so I is assumed to be equal to the observed 
width d of the jet. If along the LOS no significant changes of the quantity F fii(Te) 
occur with respect to the quantity n^, we can separate the above equation into two 
components: the emission Ije, observed at the location of the jet, and the emission 
Icor observed at the same location in the corona but before the jet. The Icor quantity 
will be given simply by 


FfniTe) 


while the I jet quantity can be written as 


r 




( 4 ) 


ijet - FfiliTg) 




'-‘F 7 


K 


dz + jl[^2^ne_cor + ngjetfdz + 


2 j 
//2 ^e_cordZ 


( 5 ) 


Because in our analysis the background pre-jet intensity was subtracted, the result¬ 
ing intensities correspond to the quantity I jet-Icor- Thus for the two filters that were 
used we can write 


jet FoidAl_Poly F^i_pgiy(Te') ' ^ ‘ (Fc jet ^ —j^t ' ^e_cor) 

(Ijet - For)Al_ _Mesh F'Al_Mesh^e) ' ^ jet J" 2 ’ ^e_jet ' f^e_cor\ 


( 6 ) 


Given the observed intensities, the only two unknown quantities in the above two 
equations are the densities Uejet and ne_cor, which can then be determined. Densities 
derived for all 18 events are shown in Figure [5bl average values of Uejet = 1.5 + 
0.1 • 10^ cm“^ and Ue cor = 1.6 + 0.2 ■ 10* cm“* were obtained. This means that 
plasma densities in the region crossed by coronal jets, Ue cor + ^ejet, are on average 
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two times higher than the plasma in the surrounding eorona Uejor 


3.4 Energy budget 


Given the extension veloeity, eleetron density, and temperature of plasma involved 
in jets, and eombining these parameters with other geometrieal and kinematieal 
parameters derived direetly by the XRT sequenees (like the jet width, length, ete.), 
it is possible to estimate different possible energies assoeiated with these events. 
As mentioned in the Introduetion, the oeeurrenee of jets is likely due to magnetie 
reeonneetions. The energy released during a jet-like e ruption may yary f rom event 
to event. Following the reeent analysis performed by Pueei et aP 1 2013 1. the total 
energy flux F provided to the eorona ean be expressed as a sum of the fluxes of 
kinetie energy potential energy Fp^,, enthalpy energy Fenth, wave energy F„ave, 
and radiative energy Frad'- 


F ~ Fkin “ 1 “ Fpot “ 1 “ Fgfitfi -|- F.^ave “ 1 “ Ff^ij 


_ 'y _ 1 

cm • 5 ]. 


(V) 


» Al_Poly/Al_Mesn 

I AI_Poly/C_Poly or ALPoly/TLPoly 
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Pucci et Ql. 2013 




In the above expression, the 
energy transfer due to ther¬ 
mal eonduetion is negleeted 
with the following rationale. 

The timeseale for eonduetion 
time is expeeted to be approx¬ 
imately two minutes, a time 
that is shorter than the av¬ 
erage jet lifetime of around 
ten minutes and mueh shorter 
than the jet radiative eool- 
ing time 12 to 40 hours [See 
Shimojo and Shibata . 2000 1. As 
a result, heat loss by thermal 
eonduetion in general should 
be eonsidered. Nevertheless, as 
mentioned before, w e eoneluded from o ur t emperature analys is (in agreement also 
with reeent results by Nistieo et al.l201 ll and Pueei et al .120 13b that the jets maintain 
a eonstant temperature along their height, so the temperature gradient is negligible, 
and thermal eonduetion along the jet ean be negleeted. 


Ng [ 1 0® cm ® ] 

Fig. 6: Plot showing possible dependenee be¬ 
tween jet density and temperature. The first 
three series of data (eirele, square, and tri¬ 
angle symbols) are those obtained here with 
different filter ratios, wh ile other data (dia - 
monds) were obtained bv IPueei et al.l 11201311 . 


The different energy fluxes have been eomputed as follows. The estimate of the 
kinetie flux Ftm is straightforward: 


kin 


1 3 
2PV 


[erg ■ cm ^ ■ s 


( 8 ) 
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where p = nefrip represents the mass density of the ejected plasma, and v is the 
outflow speed of the plasma. The potential flux Fpo, can be estimated as 

Fpot =p-gQ-l-v [erg ■ cm~^ ■ 5 “^] (9) 

where gQ is the solar gravitational acceleration (g© = 274.13 m • at the surface), 
and I is the height (maximum vertical extension of the jet). The enthalpy flux Fgnth 
depicts the thermodynamic energy of the jet and is given by 

y -2 -1 

Fenth =- 7 • P ■ V [erg -cm ■ s ] (10) 

7- 1 


where y represents the ratio of the specific heats (since the coronal plasma can be 
approximated as a monoatomic gas, we assumed y = |), and p is the thermal gas 
pressure {p = n ■ kp ■ T where kg is the Boltzmann constant). The (Alfven) wave 
energy flux F^vave represents the energy propagating outward with the jet; this is 
due to Alfven waves possibly being excited when the jet-producing reconnection 
occurs when field lines change their configuration as the y relax. Evidence fo r such 
Alfven waves appears in the jet movements observed by ICirtain et al.l KOOTII . This 
quantity can be estimated as 



_ n _ 1 

[erg - cm • 5 ] 


( 11 ) 


where B is the magnetic field strength. The quantity ^ represents the amplitude of 
non-thernial mot ions in the plasma: here we assumed ^ based on the results by 
Kim et al.l 11200711 . who found non-thermal velocities in jets ranging from 57 km s“^ 
up to 106 km s“^ . The radiative energy flux, Frad, represents the energy output due 
to radiative emission and can be estimated as 


rad 


nl ■ xT" ■ I [erg ■ cm ^ ■ s '] 


( 12 ) 


where following analytical approximations bv iPucci et al.l 11201311 . we assumed a 
0 andx = 10“^^'^^ 'N\i\\xT“ in W m^. 


The resulting averaged energy fluxes for all the 18 jets analyzed here are F^m = 
1.36 • 10^ erg cm~^s“', Fpot = 0.78 • 10^ erg cm“^s“\ and Fenth = 7.15 • 10^ 
erg cm“^s“', while the resulting F^ad component has values that are two orders 
of magnitude lower than the first three components, hence are neglected. The av¬ 
erage value of the wave energy flux is in the range Ft^ave = 0.61 • 10^ - 1.88 ■ 10^ 
erg cm“^s“', depending on the assumed values for non-thermal motion speed, and 
was computed by assuming a field strength at the base of the jet B = 2.8 G (from 
Pucci et al.l 2013t). The field in p olar regions ma y be even higher, such as 5 G 


[IPucci et al.L 1201311 or even 10 G lllto et al.l. 1201011 . The wave energy flux is there¬ 


fore probably not negligible, but its value is very uncertain due to the uncertainties 
in the magnetic field and non-thermal velocity values. As a result, in the following 
discussion we focus on the main sources of energies provided by the jets that we 
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Fig. 7. Left panel (a): The ratio between the Fg„thlPkm fluxes (j-axis), computed for the 
18 jets studied here (x-axis). Right panel (b): Plot of the enthalpy flux (on the x-axis, 
logarithmic scale) versus the kinetic flux (y-axis, logarithmic scale). Data from our analysis 
are plotted he re (black circles), together w ith t hose (red and blue circles) from the studies 
performed by Shimoio and Shib^ I 2000ll and Pucci et ah 1 2013 1. 


were able to measure independently with XRT data: the kinetic and the enthalpy 
energy fluxes. The possible energy contribution in coronal holes due to waves will 
be discussed in the Conclusions. 


4 Results 


Here we first investigate possible correlations between different physical param¬ 
eters of the analyzed jets. Figure 0 shows measured densities plotted against the 
corresponding measured temperatures for different events. This figure suggests that 
the temperature of the ejected jet plasma depends on the density of the jet for lower 
dens ity (hence fa i nter) e vents. To increase the statistics, this plot also shows results 
from Pucci et al. 1 2013 1. 


The distribution of points suggests a possible general trend: smaller, fainter events 
(those with a density lower than ~ 2 • 10^ cm“^) have a higher plasma temperature 
(higher than ~ 2 • 10® K) and vice versa. Applying a exponential decay fit of these 
points gives a Spearman rank coefficient = 0.21. This possible dependence is 
very interesting from the perspective of the estimation of total heating being pro¬ 
vided to the corona by jet events: If we consider that fainter (less dense) events 
are possibly associated with the ejection of hotter plasma, this could lead us to con¬ 
clude that very high temperature jets (T > 5 MK) are faint enough to have negligible 
emission with respect to the background corona, thus becoming barely recorded or 
not at all by XRT. It is therefore possible that the total contribution by polar jets to 
the coronal heating has been underestimated so far. 

From the estimated energy fluxes, we conclude that the largest portion of energy 
is provided by the enthalpy flux. The total amount of enthalpy-plus-kinetic energy 
fluxes has an average value of Ftot = 8.5 + 1.8 • 10® erg ■ cm~^ ■ for all the 
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events analyzed here. In partieular, the plot in Figure |7^ sho’ws that the enthalpy 
energy fluxes for all the events reported in our work are on average more than 
a faetor ~ 3 greater than the kinetie energy fluxes. As ean be seen, these results 
are also eonfirmed by data already published in the literature for other events an¬ 
alyzed with other instruments. This implies that the total magnetie energy is not 
equally eonverted into kinetie, potential, and thermal energies, as expeeted from 
the energy equipartition prineiple, usually assumed to hold for low (5 plasmas. 
Moreover, enthalpy and kinetie 
energy fluxes show a gen¬ 
eral relative eorrelation for all 
the events, as demonstrated by 
the plot in Figure [70, show¬ 
ing the enthalpy (x-axis, log- 
arithmie seale) versus the ki¬ 
netie (y-axis, logarithmie seale) 
energy fluxes derived for our 
event and those d e rived by 
Shimojo and Shibata I 2000ll to 


Fig. 8: Magnetie field strengths ealeulated via 
energy eonservation {Brec) for different jets. 

inerease the statisties. This plot 

shows a elear linear trend (in logarithmie seale): a fit to these points gives a general 


trend Fki„ = (0.26 + 0.02) • F, 


enth 


0.98 


with a Spearman rank eoeffieient = 0.92. 


The jet temperatures obtained with the ratios between different filters are on av¬ 
erage around ~ 2 MK, i n quite good agreement with previous results in the liter- 
ature : IPueei et al.l 11201311 u sed similar XRT data to find jet temperatures of ~ 1.8 


MK. ICulhane et al.l 1200711 studied temperatures and outflow speed of two jets us¬ 
ing data from t he Hinode EIS Speet rometer and derived temperatures between 0.4 
MK and 5 MK. Nistieo et al. [ 2011 1 studied jets observed in the EUVI and CORl 
eoronagraph teleseope onboard the STEREO mission, and by using the filter ra¬ 
tio teehnique between the between the EUVI images aequired with the 171A and 
195A bandpass filter s , they found values between 0.8 MK and 1.3 MK. In another 
work. lPosehek et al.l 11201011 observ ed jets with Hinode/EIS and found temperatures 
around 1.4 MK. Madjarska 1 201 ill performed a multi-instrument analysis of a jet 
using SOHO/SUMER, Hinode/XRT and EIS, as well as Stereo EUVI data and de¬ 
rived temperatures of about 12 MK for the bright point at the base of the jet and 
temperatures in the range of 0.5 MK and 2 MK for the jet plasma . All of these mea- 


rpi 

surements for the jet plasma (i.e., not ineluding iMadjarskal 11201 111 base-bright-point 
temperatures) are in general agreement, where slight differenees most likely result 
from slight variations in in speeifie properties of different jets. 


Given the values for the plasma temperature, density, and outflow veloeity of the jet, 
the magnetie field strength at the reeonneetion point ean be estimated by assuming 
that the whole enthalpy and kinetie energy eomponents are eoming from eonversion 
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of magnetic energy via the reconnection process, so that 




Brec = V^TtOoV^ + IhUbT). 


(13) 


Applying this formula to data from our 18 jets, we find an average value of Brec - 
2.3 + 0.5 G (see FigurelH). This value has to be considered as a lower limit estimate 
because our computation did not include the wave energy or the possible presence 
of turbulent kinetic energy. This average field strength corresponds, for instance, to 
the value one may expect at the altitude h ^ 5.3 x 1 0^ km (i.e., 7.3 arcsecs) for the 
field B{h) = Bo(dlh)^ [see, e.g., Aschwanden . 2013ll of a small scale dipole with its 
strength at the photos pheric surface Bn = 70 G [as observed at the base of coronal 
bright points, see e.g. iHuang et al.L 1201211 buried at the depth d = 10^ km below the 
solar surface. 


5 Summary and discussion 


A total of 18 polar jets observed by Hinode XRT were selected and analyzed to 
derive morphological, kinematical, and physical parameters. The main results of 
this work can be summarized as follows: 


• The obtained jet outflow speed, averaged over all the events, is ^ 160 + 30 
km s“^. 

• The average jet temperature was found to be 1.8 ± 0.2 MK. 

• Average electron density values of ne_cor = 1 . 6 + 0.2 • 10 ^ cm“^ and ne_cor+ne_jet = 
3.1 + 0.3 • 10 ^ cm“^ were obtained, thus plasma densities in the region crossed 
by coronal jets were found to be on average two times higher than the plasma in 
the surrounding corona. 

• The derived densities and temperatures suggest an anti-correlation between these 
two parameters, with fainter (less dense) jets being hotter and vice versa, thus 
suggesting that the total energy flux provided to the corona could be underesti¬ 
mated. 

• The largest portion of energy flux provided by jets is due to the enthalpy flux, 
while the kinetic energy flux is on average about three times smaller. The Alfven 
wave energy flux may also be substantial, but its value is uncertain. The radiative 
and conductive losses are negligible. 


Given the average energy flux for the observed jets, it is also very interesting to 
estimate the resulting total energy rate provided to the corona due to the jets and to 
compare this with the energy rate required for coronal heating. Coronal holes typi¬ 
cally have a fractional area of 0.02-0.06 of the total solar surface llSvalgaard and Kamide . 
201311 . As a result, if we assume that Ac// is the total area covered by coronal holes 


on the Sun, A jet = n ■ is the typical cross sectional area of a jet with observed 
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width d = 2Rjet, Tjet is the average jet lifetime, and that fjet is the average polar 
jet frequeney, then the total energy Eje; provided by polar jets to the eorona ean be 
estimated as 


Ejet {Fenth Fkin Fpot Fy^^ive') 


Ajet 


^CH 


'jet J jet 


fjeti 


(14) 


where the other energy fluxes have been negleeted owing to their small eontribu 
tions. In this eomputation, we assume i n partieular ^ 30 d av~^ = 3.47 • 10“^ s“' 


eorresponding to a singl e eoronal hole [ Saveheva et al.l 2007h . then t jet - 10 min 


600 s, Rjet = 8-10^ em IlSaveheva et al.L 1200711 . Ar-w - 0.04 - A, 


: 2.43 • lO^' em2, 

and {Fenth + Fun + Fpot + F„ave) = (V.15 ■ 10^ + 1.36 • 10® + 0.78 ■ 10® + 1.24 • 10®) = 
1.05 ■ 10^ erg cm“^ s“^. Here we used an average value for F^ave, based on the range 
of estimates given in Seetion 3. The resulting total energy rate provided to the 
corona averaged over the 18 events studied here turns out to be Ejet = 0.18-10'^ erg 
cm“^ s“'. This value is more than two orders of magnitude less than the energy rate 
required to h eat a polar coronal hole, whi ch is estimated to be on the order of 6 • 10® 
erg cm“^ s“' llWithbroe and Novesll 197711 . It therefore appears that the total energy 
provided by polar jets is not sufficient to heat the solar corona, unless (as mentioned 
above) we are missing a lot of high-temperature jets that are not observed by XRT 
or a very signific ant quantity of turb ulent kinetic energy, which is not included in 
our computation. Moore et"^ [ 201 ill argue, however, that if spicules are driven by 
the same processes that drive coronal jets, then, owing to the abundance of spicules, 
the sum of the energies of coronal jets and those of spicules migh t be sufficient fo r 
powering the corona. Our results also agree with recent work by Yu et al.l 1 2014 1. 
who estimated the contribution of coronal jets to the solar wind energy flux and 
derived an energetic output of ~ 1.6% compared to the total energy flux of the solar 
wind. 
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Table 1 

Appendix: Jet parameters derived with the Al Poly/Al Mesh ratio. 


Number 

Dataset 

jet width [km] 

Outflow speed [km s 

Te [MK] 

Ne [10^ cm-3] 

1 

1 

8900+750 

180+16 

1.9+0.7 

3.4+0.2 

2 

1 

5200+750 

240+33 

2.0+0.4 

4.2+0.3 

3 

1 

8900+750 

150+15 

1.7+0.2 

4.2+0.1 

4 

1 

2200+750 

120+1 

2.5+0.6 

5.6+1.1 

5 

1 

3700+750 

150+8 

1.2+0.4 

7.9+2.4 

6 

2 

5200+750 

130+5 

4.6+0.5 

0 .6+0.1 

7 

2 

3700+750 

160+14 

1.4+0.3 

4.5+0.1 

8 

2 

4400+750 

170+10 

3.3+0.4 

1 .4+0.2 

9 

2 

2200+750 

120+12 

4.0+0.6 

2.4+0.8 

10 

2 

12600+750 

270+23 

1 .6+0.1 

2.0+0.4 

11 

2 

4400+750 

130+29 

1.7+0.5 

3.8+1.1 

12 

2 

3700+750 

140+10 

1 .6+0.6 

2.5+0.5 

13 

2 

6000+750 

130+17 

1.5+0.1 

4.1+0.8 

14 

2 

8100+750 

100+13 

1 .7+0.2 

2.6+0.3 

15 

2 

8900+750 

220+4 

1.6+0.3 

3.6+0.2 

16 

2 

8900+750 

190+10 

1.9+0.4 

5.3+1.6 

17 

2 

3700+750 

140+5 

1.5+0.4 

3.4+0.6 

18 

2 

6700+750 

120+4 

1 .6+0.1 

2.9+0.9 
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